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The effect of a narrow-band isotropic stochastic GW background on pulsar timing and astrometric 
measurements is studied. Such a background appears in some theories of gravity. We show that 
the existing millisecond pulsar timing accuracy (~ 0.2 /xs) strongly constrains possible observational 
consequences of theory of massive gravity with spontaneous Lorentz braking [1], essentially ruling 
out significant contribution of massive gravitons to the local dark halo density. The present-day 
accuracy of astrometrical measurements (~ 100 /ias) sets less stringent constraints on this theory. 
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I. INTRODUCTION 

Spectacular progress in observational cosmology, espe- 
cially in measurements of CMB radiation, has challenged 
our understanding of the Universe. The standard cosmo- 
logical ACDM model based on GR is confirmed by ob- 
servations with high accuracy [2, 3]. This model requires 
the present Universe to be dominated by dark matter 
and dark energy of unknown nature, so the modification 
of gravity at large distances could provide alternative de- 
scription of the Universe. There are several theories with 
infrared modification of gravity based on quite different 
grounds (e.g. [4, 5, 6, 7, 8, 9, 10, 11, 12]). Among these 
possibilities, recently developed theories of massive grav- 
ity with violated Lorentz invariance [11, 12, 13] appear 
to be theoretically attractive and have interesting phe- 
nomenology (see the recent review [14]). In particular, 
in theory of massive gravity [13] , the Lorentz invariance is 
spontaneously broken by the condensates of scalar fields, 
which allows to avoid problems of strong coupling and 
ghosts that are unavoidable in Lorentz-invariant theories 
with massive graviton. 

Dubovsky [13] constructed a theory where gravita- 
tional waves (GWs) are massive while linearized equa- 
tions for scalar and vector metric perturbations, as well 
as spatially flat cosmological solutions, are the same as 
in GR. In this theory an extra dark-energy term appears 
in the Friedmann equations suggesting an unusual expla- 
nation to the observed accelerated expansion. In addi- 
tion, massive gravitons could be produced in the early 
Universe copiously enough to explain, in principle, all of 
the cold dark matter [1]. A distinctive feature of GWs 
produced by cold massive gravitons is a very narrow fre- 
quency range of the signal (Ai/ jv ~ 10 -6 ) as determined 
by virial motions of cold gravitons in the galactic halo. 
The central frequency itself is model-dependent, but GW 
emission from known relativistic binary systems place 
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an upper limit on the frequency v 3 x 10 5 Hz. At 
lower frequencies, the amplitude of the GW background 
could be of order h ~ lO" 10 ( 3xl0 J 5Hz ) [1] assuming the 

density of gravitons matches the conservative estimate of 
the dark matter local density pdm = 0.3 GeV cm _3 [15] . 
Clustering of GWs on ~ kpc scales constrains de Broigle 
length of massive graviton accordingly thereby placing 
a lower limit of ~ 10 _8 Hz. This leaves out a region 
~ 10~ 8 Hz < v <^ 3x 10 -5 Hz for the allowed frequency 
of GW associated with massive gravitons. The amount 
of GW signal in the frequency range ~10~ 5 — 10~ 6 Hz is 
further constrained by the tracking data for the Cassini 
spacecraft [25] (see Fig. 1). 

The aim of this note is to show that the amount of 
(almost) monochromatic GW in the entire allowed re- 
gion can be strongly constrained from the existing pulsar 
timing data and astrometric measurements, essentially 
ruling out any significant contribution due to massive 
gravitons to the density of galactic dark matter. The 
propagation of electromagnetic waves from a remote as- 
tronomical source in the presence of a GW background 
causes an excessive noise in pulsar timing [16, 17] and 
alters stochastically the apparent position of the source 
[18, 19]. So, high-precision pulsar timing and astrometry 
of distant sources (for example, quasars) can be used to 
constrain the amplitude of the possible GW background. 



II. CONSTRAINTS FROM PULSAR TIMING 

Pulsar timing was suggested in the late 1970s [16, 17] 
as a tool to detect or constrain the local GW background. 
A GW travelling through the Solar system affects the ob- 
served frequency of a pulsar resulting in anomalous resid- 
uals in the time of arrival (ToA) of pulses [20] . Because of 
unrivalled rotational stability, timing of millisecond pul- 
sars is particularly well suited for detecting GWs [21]. 
The conventional technique of the stochastic background 
measurements using pulsar timing [22] assumes correlat- 
ing ToA residuals of several pulsars. Using this method 
has yielded upper limits on the low-frequency broad-band 
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stochastic GW backgrounds [23]. 

A narrow-band GW background produces an exces- 
sive noise in pulsar timing at the corresponding fre- 
quency. The rms of timing residuals of even a sin- 
gle pulsar can put an upper limit on the GW back- 
ground amplitude in the frequency range between the 
inverse of the pulsar timing data time span T (typi- 
cally several years) and inverse time of the pulsar sig- 
nal accumulation(~ hours). The Parkes Pulsar Tim- 
ing Array (PPTA) includes several pulsars with current 
rms residuals r ^ 0.2 /xs (0.12 /xs, 0.19 /xs and 0.17/xs for 
J0437-4715, J1713+0747 and J1939+2134, respectively) 
[21]. 

In the weak field limit of Dubovsky ct al. theory, the 
equations of motion are the same as in Einstein's GR and 
we can therefore employ the results of GR calculations 
on the expected effect of a local GW on the observed 
frequency of a pulsar. For the GW power spectrum per 
logarithmic interval (as defined by Eq. (18) of [24]) we 
restrict ourselves to a (5-like function at some k = 1~kv jc: 



P h {k') 



P ,k < k' < k + Sk 
0, otherwise 



(1) 



For this power spectrum, the mass density of GWs is [24] 
Pgw = (lfor(T)- 1 c i kPo6k (2) 

providing necessary connection to the total power PoSk. 

The observed pulsar ToA rms variation r gives an up- 
per limit on ToA dispersion due to GWs and therefore 
translates into the following upper limit on P^Sk (see Ap- 
pendix A): 



P Sk < 3r 2 /fc 3 c 2 , 



and therefore 



p GW c 2 ^ (167rG) _1 3r 2 fc 4 c 4 = 3G~ 
« 2.5GeV • cm" 3 . 



VrV 



(3) 



(4) 



3 x 10" 5 Hz 



0.2 /xs 



The upper limit corresponding to r = 0.2 /xs is plotted in 
Figure 1 as a function of v and is clearly lower than what 
is needed for massive gravitons to be the dominant com- 
ponent of the local dark matter at regions unconstrained 
by the Cassini data. 



III. ASTROMETRIC CONSTRAINTS 

The astrometric effect is different for the light propa- 
gating across a region of space with enhanced density of 
massive gravitons, e.g. a dark halo of a galaxy or galaxy 
cluster (the en-route effect), and for stochastic change 
of the position of the observer immersed in the mas- 
sive graviton halo (the local effect). The former smears 
out the visible size of a distant source, while the lat- 
ter changes stochastically the angular separation between 
different sources on the sky. 




GW frequency, Hz 



FIG. 1: Astrometric (dot-dashed) and pulsar timing (solid) 
constraints on the overall energy density of a stationary 
isotropic background of monochromatic GWs as a function of 
the frequency v in the range ~ 10 _8 Hz <K~3x 10 _5 Hz 
allowed by binary pulsar GW emission and DM clustering 
constraints (see Introduction). The thick dashed line corre- 
sponds to the local DM energy density of 0.3 GeVcm -3 . The 
constraint in the lower right corner of the graph is set by the 
Doppler tracking of the Cassini spacecraft [25] . 



In astrophysically relevant cases the en-route effect is 
too small to be detected at the present level of accuracy of 
astrometric measurements. A very generous upper limit 
on the stochastically fluctuating change in the observed 
position of a distant source may be estimated as aq, ~ 
/i max $L, where is the angular size of the halo on the 
line of sight and h max is the maximum amplitude of GWs 
comprising the halo (which, due to the Gaussian nature 
of these fluctuations, is essentially the same as the rms 
amplitude h c = \fP~h that can be estimated from the dark 
matter density). 

Contrary to naive expectations, in GR the light ray 
deflection does not execute a random walk and does not 
show the ~ \/N growth of the deflection. Instead, for 
traceless tensor perturbations travelling with the speed 
of light, only the gravitational wave field at emission and 
detection points matter [18, 26, 27]. The relative change 
Aty /iff in the angular separation between two sources due 
to the local effect is also of order of the GW background 
amplitude h c . At h c ~ 10 -11 — 10~ 10 , as Dubovsky et al. 
model suggests [1], this would yield ~ /xas jitter in the 
angular separation for a couple of sources across the sky. 
Such jitter can be discovered in the future astrometric 
space experiments like SIM [28]. 

Present-day astrometric accuracy a* can be estimated 
as that of the radio VLBI-based ICRF (International 
Cosmic Reference Frame) [29, 30], which involves more 
than 200 reference radio sources determining the celes- 
tial coordinate frame. The ICRF sources are observed 
for many years, and the accuracy of determination of 
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source coordinates on the sky relative to this frame may 
be used as a measure of the angular separation stability. 
The best present-day accuracy of 100 /zas (at 1 a level) 
[31] means A*^5x 10~ 10 . 

For the adopted power spectrum (1), this accuracy 
translates into the following upper limit on the GW mass 
density (see Appendix B): 



3ttz/ 2 ct| 
4.2 x 10 3 GcV ■ cm" 3 



(5) 



3 x 10- 5 Hz 



0\p 



100 /ias 



Other limits for this value as a function of frequency v 
are summarized on Fig. 1. 



IV. CONCLUSIONS 

We have shown that the existing data on the millisec- 
ond pulsar timing stability set a tight upper limit on the 
narrow-band GW background amplitude at frequencies 
v ^ 10~ 5 Hz. This limit can be used to severely bound 
the amount of massive cold gravitons which can poten- 
tially produce a strong narrow-band GW background [1] . 

The present-day astrometric constraints are less re- 
strictive than the timing ones at considered frequencies. 
However, both are still far above the tightest constraint 
set by the Doppler tracking of Solar system spacecrafts 
in the frequency range v > 10~ 6 Hz [25]. 



of e, k' and choice of the polarization gauge; in the above 
formula, summation of circulary polarized modes s = 1, 2 
is implied and 'c.c.' stands for the complex conjugate. 

The time of arrival residual R(t) is the integral of rel- 
ative fluctuation Af(t)/f w.r.t. time: 



(A2) 



R(t) = / dtAf(t)/f 



n 

and its autocorrelation function is 

T 



R 2 {t) = lira T- 1 J dtR{t)R{t + r). (A3) 



In particular, plugging the above expressions into each 
other and using statistical properties of h s (k, t) as set by 
Eq. (18) of [24] one obtains for the ToA residual rms 



r A = R 2 (0) = (3c 2 )" 1 / dk' P h (k')k'-' 



(A4) 



Applied to the power spectrum (1), this imediately gives 
P Q Sk = 3r 2 k 3 c 2 (A5) 
and corresponding limit (4) on the GW mass density. 



APPENDIX B: ASTROMETRIC FLUCTUATIONS 
IN NARROW-BAND GW BACKGROUND 
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APPENDIX A: PULSAR TIMING IN 
NARROW-BAND GW BACKGROUND 

To calculate the effect of a narrow-band GW back- 
ground on pulsar timing one can use the formalism pre- 
sented in Baskaran et al. [24] starting with their ex- 
pression (9), which describes the effect of a GW with 
wavevector k' and polarization tensor py on the observed 
rotational frequency / of a pulsar at position e on the sky. 
Frequency fluctuation due to modes with all s and k' is 

Af(t)/f = J dk' [h s (k',t)g s ( PlJ1 k',e)+c.c] , (Al) 

where g(pij, e, k') is given by (9) of [24] (setting the pa- 
rameter e — > 0) and depends on the relative orientation 



To calculate the fluctuation An in the position of a 
source n we write the astrometric effect due to a sin- 
gle GW mode in the formalism of [24] similarly to Ap- 
pendix A: 



An(i) = f s (k,n)/i s (k,t) + c.c. 



(Bl) 



where conventions following (Al) are respected and fac- 
tors i s (k, n) can be found by appropriately rotating the 
results of Pyne et al. [32] : 



f 1,a (k,A) 



n. 



k,n 



±i 



n, k 



}■ 



(B2) 



As An _L n, the former has only two independent com- 
ponents. For a pair of sources n, n' it is natural to choose 
these components to be along the great circle connecting 
the pair and perpendicular to it. Using the above formu- 
lae and the statistical properties of h s (k, t) as introduced 
in Eq. (18) of [24] one obtains for the correlation matrix 



Anil (n,i) An" (n'.tf) Anll(n,t)An- L (n , ,i')\ fml 
An x (rM)A7iH(n',t') An- L (n, t)An 1 (n', if)) ~ ( } 



l A o\ 2 * D , , , l (l 

- I r. -i cos — Pat) + - „ 
6 V0 1) 2 cW 6 V° 



o A . 2 * 



sin 2 -P C (2T) 
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where * = Z(n, n'), T = t' - t, 2T = t' +t and P c (i) is 
the cosine Fourier transform of the GW power spectrum 
P h (as defined by Eq. (18) of [24]): 

P c (t) = f dkk' 1 P h (k) cos ckt; (B4) 

for a (5-like power spectrum (1), P c {t) = k~ 1 P( j 8k cos ckt. 

The autocorrelation function of the principal observ- 
able, the fluctuation is 

(A* (t) A* (f )) = I sin 2 1 [P c (r) - P C (2T)] . (B5) 

Independent observations average away the second term 
in the square brackets above in accord with the station- 
arity of the problem. In particular, the rms value (r = 0) 
in this case is 



(A *»> = Isi4p,,(») = — s,»--; 
hence the estimate (5). 



P 6k 



(B6) 



The gravitational wave field is Gaussian and so is its 
linear transform A\I>(i) (cf. Bl). Therefore (B3) readily 
gives the respective distribution functions. The results 
generalize naturally to multiple observables. For a set 
of m angular distances ^i(U) = + A\tj(ij) between 
source pairs (kj,n.i), the probability density y>(A\P) of 
the observable vector 

A* = (A* 1 (t 1 ),...,A* m _ 1 (t m _ 1 ),A* m (t m )) T 



is 



<^(A*) 



exp 



iA^M 1 (ni,...,!^; 



A* 



v /(27r) m dct M 
with elements of M given by (<E>y = /(n, + k i? hj + Ikj)): 



1 *j , 
My = - sin — sin -± cos §ijP c {U -tj). 
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